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Abstract 




The experimentally observed ion production rate in a cesium thermionic converter with 
tungsten electrodes is 3 to 4 times the rate predicted by the Saha-Langmuir equation. The ion 
emission rate is found to be dependent on field and corresponds to an exp (e-yel/kT) (Schottky) 
relation. The Schottky effect is unambiguous only at very small spacings (3 to lip) and at 
conditions of low ionization probabilities. Agreement between ion current and simple space- 
charge theory is observed. This agreement and the departure from the Saha-Langmuir relation 
suggest that the ions are formed at localities of higher work function than the electrons. 

The difference in work function for ion and electron production thus reduces the effectiveness 
of ions in reducing the space charge of electrons at the electron- saturation condition. Obser- 
vations on the coupling of ion motion to particle arrival rate are also made. 


Introduction 




The "ideal" passive-mode thermionic converter (one in which volume ionization does not 
occur) is still an important objective in thermionic conversion. It is probable that surface 
ionization will play a significant role in thermionic conversion if an effective passive-mode 
operation is to be achieved. This paper presents observations on ion production probability, 
the Schottky effect on ion generation, space-charge effects, and the coupling of ion motion 
with arrival rates. 


Symbols 


A 24, constant in eq. (A6) 

B 1.5, constant in eq. (A6) 

E field 

e electronic charge 

3F function 

h Planck 1 s constant 

J current 

J Q zero field satuation current 

k,k* Boltzmann constant 

m mass 

n number density 

n a number density of neutral atoms 

He number density of electrons 

n+ number density of positive ions 

p pressure 

Q partition function 

Qa partition function for atoms 

Qe partition function for electrons 

Q+ partition function for positive ions 

q electronic charge 

r a atom reflection coefficient 

r + ion reflection coefficient 

T temperature 

T cs cesium reservoir temperature 

emitter temperature 


€ 0.577, constant in eq. (A6) 

0 angle 

K 2.492, constant in eq. (A6) 

m free energy 

^ free energy of atoms 

Pe free energy of electrons 

lip plasma potential 

free energy of ions 
v particle flux 

v a atom flux 

v e electron flux 

v e p electron flux from plasma 

v es electron flux from surface 

v + ion flux 

p dimensionless density 

a cross section 

o m total cross section 

ff +,m energy-dpendent momentum, transfer 
cross section 
9 work function 

9^ emitter work function 

92 collector work function 

X dimensionless distance 

i|r dimensionless potential 


V voltage 

V i ionization potential 

x distance 

a ratio of ion particle flux to atom particle flux 
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Analytical Relations 

The relations given in Appendix A relate to ion emission and ion transport and are com- 
monly used in thermionics. To facilitate the review of the experimental results, however, a 
"brief resume of the Saha-Langmuir equation is presented. 

The Saha-Langmuir equation is based on the equilibrium between an ionized gas (plasma) 
and the containing surface. The flux of the electron gas from the metal containing surface 
must be equal to the electron flux from the plasma. For the sake of simplicity, it is assumed 
that the body of the plasma and the plasma near the surface are at the same potential as the 
surface (pp = ecp) . The Saha equilibrium equation is then used to define the relative particle 
concentrations in the plasma. The expression for the free energy of this system is 

p a = P+ + Pe + eVi (l) 


Substituting equation (A9) in equation (l) results in 

n+n e Q+Qe ( 2 it meklY 5 / 2 / eV i\ 

n a " Qa ( h2 j SXP \ kT ) 

In terms of particle flux, since v = n^/iS'L., 

™2:rtm 



( 2 ) 

(3) 

(4) 


The equivalent development of the Richards on- Du shman equation for the electron gas in the 
metal with its surroundings^ results in 


4it m e (kT) £: 



( 5 ) 


Equating the surface and plasma electron fluxes (v es = v ep ) and solving for the ratio of 
ion to atom flux yield the Saha-Langmuir equation 


v a 


■VkT/2jtm e ( 2 jtm gkT/h 2 ) 


3/2 


exp (- eV^/kT) 


4inn e [(kT) 2 /h 3 ] exp - (ecp/kT) 


( 6 ) 


[e(q> - Vi)] 

= 1/2 exp[— J 

At equilibrium, the flow to the surface and away from the surface for each particle is, 
of course, equal. Many refinements to the Saha-Langmuir equation exist. For example, con- 
sideration of particle reflection (treated by Copely and Phipps 2 ) results in 


v + (1 - r+) x fe(<P - Vi)] 

^ = (1 - ra) 2 ex P [— kT J ( 8 ) 

The effect of field on surface ionization has also been treated in a survey article by 
Zanberg and Ionov, 3 in which Morgulis (1934) is credited with first observing the field effect 
on surface ionization, and Dobretsov (1936) is credited with relating the ion production 
enhancement to the change in image field (analogous to the Schottky effect). 

The simple equilibrium approach used in equation (6) can be continued to include these 
field effects by first assuming that the plasma potential is not equal to the surface potential. 
The number density of the particles in the plasma is still described by equation (2). 
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Matching the plasma fluxes to surface fluxes at equilibrium requires that the Boltzmann 
relation (eq. (A3)) be maintained. Thus the potential corresponding to the increased charge 
density. Up < ecp, or reduced charge density, pp > ecp, must match at the surface where pp 
equals eqp e ff • Since equilibrium constraints are still employed, the relative flux rates are 
yet described by equation 4. Xt follows that a plasma sheath must exist to permit matching 
the plasma boundary conditions with the surface boundary conditions. (These relations have 
been observed and evaluated numerically by Nottingham 4 in his isothermal diode paper. ) 

Since equation (4) applies, equation (6) must apply. But the application of equation (6) 
to the case where pip is less than ecp produces a strong electron accelerating field away fro® 
the surface. The work function cp is modified for the effect of the field (the Schottky 
effect - Nottingham 4 ), Equation (6) then becomes 


/kT ( Z * m e kT \ 3 /2 

I± IJEV \ h 2 / exp - (eVj./kT) 

Va ■linnet (kT) 2 /h 3 ] (exp - <pe/kT)[exp (e-^/eE/kT)] ^ 

Cancelling terms as before results in 

— = i exp [(ecp - eVi - e-y/eE) /kT] (10) 

Va, £ 

for fields that accelerate surface electrons. By analogy, a similar image force relation 
should occur for a field that accelerates the positively- charged ions. The effective work 
function barrier that determines charge flux would then be increased by an amount 

Of = e-^/eE At (11) 

and equation (10) vould become 


= - exp [(e(p - eVi + eVeE)A T ] 
v a c 


for a field that accelerates ions. 

The Saha-Langmuir equation can be extended to include "patchy" (variable work function) 
surfaces. The various formalisms, while not complex, are bulky, so they will not be repeated 
here (see Zandberg and Ionov 3 or Kaminsky 5 ). Zandberg has indicated that field as high as 

10 6 v/a may be necessary to eliminate "interpatch " effects in field-enhanced surface- 
ionization studies. 

The equilibrium-based equations appear to be reasonably supported by experiment 
(Kaminsky 5 and Zandberg and Ionov 3 ). Although considerable contradictions exist in past lit- 
erature on the field effect on surface ionization, a more recent study by Zandberg and Ionov 6 
support the e-y/eE /kT dependence. Less well established are surface -ionization studies at 
conditions encountered in thermionic conversion, that is, high- fractional cesium coverage and 
low ionization probability. 


Apparatus and Procedure 

The various features of the planar variable -spacing diode are essentially as described in 
detail by Breitwieser 7 . The 1. 5-centimeter -diameter emitter was formed from tungsten oriented 
to expose the 110 (Miller index) plane. Several crystals composed the face, but the maximum 
misorientation was only 1^ . A post-mortem examination disclosed some surface reorientation 
to the crystal faces. The body of emitter subsurface retained a 110 orientation. It is not 
clear whether the surface reorientation existed during the tests or was a result of damage 
incurred as the tests were terminated by a leak. 
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A refinement to the operating procedure previously described^ was the method of use of the 
collector guard and an associated high-gain differentially coupled amplifier. Improved current 
resolution in the microampere range resulted. 

Presentation of Data 


Figure 1 illustrates the comparison of experimental ion currents with ion currents cal- 
culated by equation (12). The evaluation of work function was determined through the 
Richardson Dushman equation (Al). The experimental data were extrapolated back to the zero- 
field condition through the use of the Schottky relation (A2). The observed data were obtained 
at conditions corresponding to field strengths up to 3.5xl0 6 V/m. The data shown were obtained 
at spacings estimated to be 3-11 p. Knudsen flow was assumed in calculating the arrival rates 
in the analytical estimates. 

Figure 2 is a typical Schottky plot taken at an emitter temperature of 1530° K, at a 
cesium reservoir temperature of 470° K, and at three interelectrode spacings. The condition 
of operation would correspond to an ion-rich condition at the emitter on the basis of the 
Saha-Langmuir equation (A4) in which the emitter work function is determined from the saturated 
current of electrons. The quantity V - 1 in the abscissa adjusts the applied voltage V to 
the surface potential difference (cp^ - cp 2 ** l)- The smallest spacing is estimated to be 3 to 
11 p on the basis of external measurement. The Schottky slope indicates a spacing of 10 p. 

The vertical lines labeled zero-field potentials were calculated by using the Langmuir space - 
charge relation for singly charged cesium ions (eq. (A8)). No electron neutralization was 
assumed. The extrapolated zero -fie Id ion current based on the Schottky analysis for the 
closest spacing data was used in the dimensionless distance relation for all spacings 
(eq. (A8) ) . The resulting values of voltage to produce the zero field are plotted as vertical 
lines. The dashed vertical line indicates the region where the electron current would be 
expected to have a significant effect on the net current, as estimated from close-spacing data. 

Figure 3 is similar to figure 2. The temperatures shown correspond to a slightly 
electron-rich condition at the emitter on the basis of the Saha-Langmuir theory and electron 
currents but slightly ion-rich on the basis of the observed currents. Again, the small-spacing 
Schottky slope (10 p) corresponds to the measured values for the closest spacing (3 to 11 p). 
The potentials required to produce a zero field at the emitter (single -charge theory) intercept 
the actual data to indicate a consistency in zero field currents. 

Figure 4 is similar to the two previous figures. The emitter temperature is 1600° K, and 
the cesium temperature is 530° K. Thus, the probability of ion-neutral collision is greater. 
Included in figure 4 are runs that represent the level of data scatter involved in successive 
(days or weeks) runs. The results are similar to results presented in figures 2 and 3 except 
that the zero field potentials do not intercept at the same value of ion-current density. 

Figure 5 is a plot of the logarithm of current as a function of applied voltage. The 
interelectrode spacing was 3 to 11 p. The emitter temperature of 2007° K and the cesium 
reservoir temperature of 309° K were selected to ensure nearly complete ionization of the 
arriving particles (exp [e(cp - Vi)/kT] > > l). It was anticipated that these conditions would 
yield information on particle arrival rates and thus clarify some aspects of the deviations of 
ion current from the Saha-Langmuir theory. The results indicated instead that a coupling 
exists between ion motion and arrival rate. 

Two theoretical curves were superimposed on the experimental data. The higher values of 
theoretical ion current were based on free -molecule (Knudsen) flow. It was assumed the flux 
rate above the cesium reservoir, as calculated by the vapor pressure relation of Heimel^ 

(eq. (A7) ) , was identical to the flux in the interelectrode space (Knudsen flow). The theo- g 
retical curve was corrected for ion space charge by the method of Nottingham and Breitwieser y 
The lower theoretical values of ion current assume that the interelectrode space is at the 
same pressure as the surrounding gas. It was further assumed that the momentum exchange was 
primarily between the surrounding gas (assumed at T cs ) and the ions in the inte relectr ode 
space. Thus the flux in the interelectrode space was reduced by the ratio of V T cs/^l* 

Qualitative agreement between the ion -current and constant -pres sure arrival rates exist at the 
low retarding voltages customarily used in thermionic converter studies (-3 to -8 V) , but at 
high retarding voltages (-30 to -40 V) the Knudsen flow provides a more adequate comparison. 
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Figures 6 and 7 are extensions of figure 5, in which the same comparisons are made at 
larger spacings. In some cases the complete development of the theoretical curve was termi- 
nated at the condition of zero field at the emitter, thus avoiding the tedious space -charge 
calculations. It should be noted the cesium reservoir temperature varied from 325° to 329° K 
in these tests; thus the slight difference in the value of the saturated ion current. All sets 
of data in figures 5, 6 and 7 show the same general trend. The observed ion currents at -30 to 
-40 volts correspond to the Knudsen predictions. Currents measured at 1 to 2 volts more nega- 
tive than the theoretical zero-field condition at the emitter correspond to the constant- 
pressure predictions. 


Discussion of Results 


Although the magnitude of the experimentally observed ion current is greater than that 
predicted by the Saha -Langmuir equation, the close agreement to the functional form of the 
equation is indeed surprising in view of the equation’s equilibrium origin. The departure 
from equilibrium, while not great for electron emission, would appear to be a major perturba- 
tion in the case of ion emission, particularly at high fields. The experimental values are 
consistently only 3 to 4 times that predicted by the Saha -Langmuir relation. Possible reasons 
for the departure are many. Sxperimental errors in temperature measurements, such as emitter 
and cesium reservoir, were explored but could not adequately explain the deviation. Arrival- 
rate discrepancies, based on the various cesium vapor pressure relations or distortion in 
local fluxes, were insufficient or in the wrong direction to explain the higher ion-production 
probability, as indicated in the particle arrival -rate studies illustrated in figures 5, 6 
and 7. 

One self-consistent possibility exists: the work function effective in producing ions 

may be higher than that related to the emission of electrons. For, example, a difference of 
1.2 e/kT (the eV equivalence of emitter temperature) would account for the difference be- 
tween experiment and theory shown in figure 1 if equal areas for ion and electron emission are 
assumed. This difference in surface potential level then could reduce the electron current 
(and thus electron density) at the zero-field condition for ion emission sufficiently to per- 
mit the use of the single-species space-charge theory. The justification for the use of this 
space-charge analysis was in part established by examining the analytical mixed-charge space- 
charge results of Goldstein-l*0,ll. The results indicate that if the ion emission density is 
about 10 times the electron emission density, the single-species space-charge theory is ade- 
quate to determine the zero-field emission condition. At low values cf X (eq. (AS)) the 
approximation improves. 

Since the results in figures 2 and 4 correspond to an ion-rich emission density condition 
(based on the application of the Saha-Langmu ir relation to the work function calculated from 
electron current), an additional barrier of 1.2 e/kT should then allow the experimental results 
to approach the required ion to the electron emission density ratio of 10. However, the condi- 
tions of figure 3 correspond to an electron-rich condition. Thus a difference slightly greater 
than 1.2 e/kT may be required to satisfy the single-species space-charge criteria. At this 
time there does not appear to be an unambiguous method of estimating the various work functions 
and area distribution suited to correlate the electron and ion emission. 

The close agreement between the zero-field ion currents at various spacing and tempera- 
tures lends credence to the space-charge analytical methods. Agreement of the data heretofore 
attained only at the close spacing (3 to Up) is now attained at a variety of spacings and 
again supports the observation of a consistent deviation from the Saha -Langmuir relation. 

Slopes in the space-charge region are obviously unrelated to the Schottky analysis. 

Previous observations by Nottingham*^ , indicating agreement between isothermal diode 
theory and ion current measurements at conditions of low ionization probability, are not veri- 
fied in this study. It should be noted that some of the present results include applied 
fields that satisfy the Zandberg- Ionov 3 criteria for nonanomalous Schottky effects in the 
presence of patchy surfaces. These data are the close-spacing data, for which all the spacing 
calculated from Schottky slopes are within the mechanically measured spacing estimates 
(figs. 2, 3, and 4) . 

Houston 13 has suggested that the reduction in ion current with spacing, observed In 
figure 3, is due to ion-mobility effects rather than space-charge effects. His analysis is in 
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part based on isothermal diode analysis and the usual continuum-transport relations. In order 
to judge the applicability of mobility analysis the mean free path must be estimated. From 
the cross-section equation (eq. (A6)), the mean free path for an ion-neutral encounter (includ- 
ing charge exchange and small angle scattering based on momentum exchange) is about 50p at 
T cs of 470° K, l£p at a T cs of 500° K, and 7p at a T cs of 530° K. (The relative velocities of 
atom to ion are assumed to be equivalent to about 0.1 eV in these estimates.) Because of the 
high accelerating potentials applied (up to 30 V) , the velocity gained from the field in one 
mean free path is greater than the thermal velocities of the particles. It therefore appears 
that the necessary criterion for conventional mobility analysis does not hold. Indeed the 
actual mean free path should be larger than the numbers previously cited since the relative 
velocity may be considerably greater than initially assumed (see kT effect in eq. (A6)). Yet it 
is obvious that collisional -scattering effects should influence ion transport. The effect is 
clearly demonstrated in figure 4, where the mean free path is the smallest, and the zero-field 
space-charge solution for the larger spacings does not agree with the close-spacing values. 

No adequate analysis as yet exists for ion transport, which includes electron flow, 
charge exchange, and velocity-dependent cross section, in the case of few-collisions. Intu- 
itively one can argue that collisional activity can be partially compensated for by electron 
space charge. To date this appears as the only other explanation of the agreement of single- 
species space -charge theory with experiment, other than the fact that ions and electrons were 
not formed at the same surface potentials at the dynamic conditions of the tests. 

The analysis of the data of the effective arrival rates in the interelectrode space, 
presented in figures 5, 6, and 7, is complex. Three particle groups are interacting. The 
ion group is at temperature Tp and moves at a velocity dictated by the temperature of the 
particles and by the difference in surface potentials. The atoms from the collector exist at 
temperature T 2 , while the atoms around the gap are at T cs . The analytical solution of the 
problem entails all the complexity of a radiation-cavity problem complicated by collisions 
that are velocity dependent and coupled to applied potentials. 

It is expedient, therefore, to look at the extremes. The relative -momentum cross sec- 
tions of the atom-atom to atom-ion collisions are about 1 to 3 (see eqs. (A5) and (A6)). Thus 
at low voltages (near zero -surface field) , the ion -neutral collision in the gap would appear 
to dominate because of the high ion density and relatively large cross section. An equilib- 
rium concentration for the cesium particles is probably dictated by momentum exchange between 
incoming atoms and ions in the interelectrode space. As discussed previously, if flux based 
on constant pressure between the gap and its surroundings is assumed, approximate agreement 
between experimental and theoretical ion arrival rates is noted at low fields. 

As the voltage is increased, the average residence time of the ion is decreased. It 
follows that the ion concentration in the gap relative to the atom concentration decreases and 
the flux balance in and out of the gap shifts toward the atoms. Knudsen flow appears to pre- 
vail. Within the framework of the questionable model assumed, the vapor -pres sure relations 
used seem to match the experimental data. 

The coupling of arrival rate with applied voltages is not belabored because of any par- 
ticular practical implications to thermionic -diode performance or as absolute verification of 
vapor pressure data but only to prevent possible errors of interpretation in some aspects of 
thermionic -converter diagnostics. 


Conclusions 


1. Surface ionization of cesium on cesium-coated tungsten surfaces is consistently 3 to 
4 times greater than that predicted by the Saha -Langmuir theory. 

2. A field-dependent ionization probability is observed following the Schottky mirror - 
image force relation. 

3. The magnitude of the ion current and the applicability of the single-species space- 
charge theory suggest that ions are formed at a surface of higher work function than are the 
electrons . 

4. At conditions of high-ionization probability the ions couple with the atoms to alter 
arrival rates on diode surfaces. 
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Appendix A 
Analytical Relations 
Richardson Dushman equation 

J 0 = 120.1 T 2 exp(- etp/kT) 

Schottky equation 

J = J eg>(eVS/KP) 


Boltzmann relation 

J = J 0 exp( -eV/kT) 

Saha- Langmuir equation 

a = — exp[(ecp - eVi + e VSs)/ kT l 

v a c 

Neutral- neutral cross section from Sheldon and Manista 14 

o = f 2ita(0)(l - cos e)(sin 0)dfi 

J 0 


(Al) 

(A2) 

(A3) 

(A4) 

(A5a) 


Neutral- neutral cross section with c(e) values of reference 14 


470 Pr at 400° K 


(A5h) 


Ion-neutral cross section from Sheldon 15 

°+,n = 2A 2 - 4AB (| - € + In k'l) + E 2 ]^ - e)ln k'T 

+ 2(ln k'T) 2 + k] +(|) 2 


e 2 a 


c, *d - ** “7 

o + n * 1400 A 2 at k'T = 0.1 eV 


\2 


(A6) 


Cesium vapor pressure relations from Heimel 3 
log 10 p (nevtons/sq. m) = - - 


0.915282 log 10 T + 12.05025 


(A7 ) 


The La ngmu ir space-charge relation^ X expanded in terms of usual variables for a 

singl y charged cesium ion 


(1. 442X10 )J 1 / 2 X T -3 / 4 


(A8) 


The symbol X as a function of + is tabulated in several references (e.g., Notti ngham 1 7 ) , 
Free energy per particle (Boltzmann statistics) 1 


p = kT In 


( jl\ 5/z a 

\2jnnkT/ Q 


(A9) 
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EXP. SAT. ION CURRENT DENSITY, A/SQ CM 



FIGURE 1. - COMPARISON OF ANALYTICAL AND EXPERIMENTAL 
ZERO-FIELD ION CURRENTS IN A CESIUM-TUNGSTEN DIODE; 
EMITTER TEMPERATURES, 135(fTO 1650° K; INTERELECTRODE 
SPACING, 3TOU|i. 



FIGURE 2. - SCHOTTKY PLOTS OF ION CURRENT 
SHOWING ZERO-FIELD INTERCEPTS BASED ON 
SINGLE-SPECIES SPACE-CHARGE ANALYSIS; 
EMITTER TEMPERATURE, 1530° K; CESIUM RESER- 
VOIR TEMPERATURE, 470° K. 



FIGURE 3. - SCHOTTKY PLOTS OF ION CURRENT SHOWING 
ZERO-FIELD INTERCEPTS BASED ON SINGLE-SPECIES 
SPACE-CHARGE ANALYSIS AT NOMINALLY ELECTRON- 
RICH EMISSION CONDITION; EMITTER TEMPERATURE, 
1500° K, CESIUM-RESERVOIR TEMPERATURE, 500° K. 




eh/V-1 /kT CS -37554 


FIGURE 4. - SCHOTTKY PLOTS OF ION CURRENT 
SHOWING ZERO-FIELD INTERCEPTS BASED ON 
SINGLE-SPECIES SPACE-CHARGE ANALYSIS; 
CESIUM PRESSURE CORRESPONDING TO 
CONDITIONS OF SIGNIFICANT ATOM-ION 
COLLISIONS; EMITTER TEMPERATURES, 1600° 
K; CESIUM TEMPERATURE, 530° K. 


VOLTAGE, V 


-30 -25 -20 -15 -10 -5 0 



FIGURE 5. - EFFECT OF APPLIED VOLTAGE ON PARTICLE ARRIVAL RATES; EMITTER 
TEMPERATURE, 2007° K; EMITTER WORK FUNCTION, 4.92 V; SPACING, 3 TO 11 |i ; 
CESIUM-RESERVOIR TEMPERATURE, 329° K; EMITTER AREA, 1.765 SQ. CM. 
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VOLTAGE, V 



RGURE 6. - EFFECT OF APPLIED VOLTAGE ON PARTICLE ARRIVAL RATES; EMITTER 
TEMPERATURE, 2003° K; EMITTER WORK FUNCTION, 4.92; SPACING, 113 M; 
CESIUM RESERVOIR TEMPERATURE, 329° K; EMITTER AREA, 1.765 SQ. CM. 


VOLTAGE, V 



FIGURE 7. ■ EFFECT OF APPUED VOLTAGE ON PARTICLE ARRIVAL RATES; INTER- 
ELECTRODE SPACING, 227 AND 452 U; EMITTER AREA, 1.765 SQ. CM. 


